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Within the model of elastic impurity scattering, we study how changes in the energy dependence 
of the electronic density of states (EDOS) N[e) around the Fermi energy ep are reflected in the 
frequency-dependent optical conductivity o{uj). While conserving the total number of states in 
A'^(e) we compute the induced changes in criuj) as a function of ui and in the corresponding optical 
scattering rate l/rop(aj). These quantities mirror some aspects of the EDOS changes but the rela- 
tionship is not direct. Conservation of optical oscillator strength is found not to hold, and there is 
no sum rule on the optical scattering rate although one does hold for the quasiparticle scattering. 
Temperature as well as increases in impurity scattering lead to additional changes in optical prop- 
erties not seen in the constant EDOS case. These effects have their origin in an averaging of the 
EDOS around the Fermi energy ep on an energy scale set by the impurity scattering. 

PACS numbers: 78.20.Bh, 78.20.-e 



I. INTRODUCTION 

Measurements of the infrared conductivity cr{ijj) as a 
function of energy uj, continue to give valuable informa- 
tion on charge dynamics in a wide range of metallic svs.- 
tems including the high- Tc superconducting cupratesBH. 
These materials have received a lot of recent attention 
because they represent strongly correlated systems which 
exhibit new physics, beyond the usual Fermi liquid (FL) 
description of electric structure. It was recognized and 
emphasized as crucial very early on, that the normal state 
properties of the cruppates are anomalous. A Marginal 
Fermi liquid (MFL)|3'QEI phenomenology was developed 
which could describe remarkably well many of the ob- 
served deviation from FL behavior of the normal state. 
An essential feature of the MFL is that quasiparticle 
weight in the single particle charge carrier spectral den- 
sity denoted by goes to zero logarithmically as the 
Fermi energy is approached. In this limit, there are 
no well defined quasiparticle poles, and the entire spec- 
tral density consists of an incoherent background which 
is due to the interactions. It is the delta function-like 
quasiparticle contribution (broadened by the interaction) 
which leads Xq, a Drudc like contribution in the optical 
conductivityQi3. The incoherent background is responsi- 
ble for the Holstein tails due, for example, to phonon as- 
sisted absorption in the well studied case of the electron- 
phonon interaction. The incoherent contribution to the 
optical conductivity gives additional information on cor- 
relation effects complementary to the Drude response. 
Both contributions are described microscopically by the 
electron self energy S(ci^) vs uj, which is the fundamen- 
tal quantity about which we would like information from 
measurements on the optical conductivity cr{ijj). For ex- 
ample, the real part (Si) of S deals with mass renor- 
malization of the quasiparticles (qps) and the imaginary 
part (S2) is related to their lifetimes. As we have just 
described, S can also lead to an incoherent background. 

One of the most striking manifestation of correlation 
effects in the high- Tc superconducting cuprates are the 



pseudogap features observed in their normal state. They 
are particularly prominent in underdoped systcrps, but 
are also known to be present at optimum dopingj. The 
precise origin of the pseudogap is not yet known and this 
remains a controversial area. Nevertheless, the experi- 
mental situation is reasonably well characterized and has 
been reviewed by Timusk and Statta. 

The pseudogap has been identified as a distinctive 
and sometimes even abrupt change in the tenuperature 
variation of_the nuclear spin lattitje-. relaxations, of the 
Kniffih|l|-S|ljiiftLl, of the dc resistivityll30, and of the specific 
heatPrr'ir'i, lifl-the frequency dependence of the infrared 
conductivityollj and tho, current voltage characteristics 
of a tunneling junctional aa-well as in angular resolved 
photo emission (ARPES)Bt3'El. This last experimental 
technique is particularly powerful and has revealed that 
the pseudogap is not constant around the Fermi surface. 
Rather it has a d-wave nature which is the same symme- 
try as is exhibited by the superconducting gap below 
in the cuprates. 

That the pseudogap has its origin in correlation ef- 
fects is not in doubt. Rather, the_issim.ia_hDBi.i±_is to be 
simply, yet accurately, describecES-ElBBOy-El. Many 
theoretical suggestions have been made. One widely held 
view is based on the so-called preformed pair model in 
which it is envisioned that the Cooper pairs exist above 
Tc up to a higher pseudogap temperature T*, but with- 
out phase coherence. The phase coherence between the 
pairs which is essential for suiiep;onductivity, sets in only 
at lower temperature T < TcCJO. In another model, dif- 
ferent but related, finite momentum pairSjare believed to 
be responsible for the pseudogap featureso. A very dif- 
fcrcnt_recent proposal is the suggestion of Chakravarty 
et aZ.EEl of D-density wave formation with attended or- 
bital currents which double the crystallographic unit cell. 
There are also pronosals encoded in the ideas of spin- 
charge separatiorOO and the pioneering suggestion of 
Andersoric3. 

The true nature of the changes that are brought 
about in the energy dependent electronic density of states 
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(EDOS) N{e) by the formation of the pseudogap remains 
unknown other than that the EDOS is depressed in some 
way. Consequently we will not address this specific case 
directly here although it is a motivating force for what we 
have done. Instead we will be concerned with a related 
but less specific issue, namely, the general question of 
how changes in N(e) around the Fermi energy will mani- 
fest themselves in corresponding changes in the frequency 
dependence of the optical conductivity <t(w) vs uj. After 
all, from an experimental point of view, it is important 
to understand what qualitative signature is to be looked 
for which correspond to microscopic changes in N{e). 

To remain as simple as possible, we will examine in 
this paper in some detail mainly a simple model for N{e) 
which consists of a constant background Ni, modified by 
two Lorentzian forms, both chosen to be symmetric about 
the Fermi energy. This assumption allows us to take ad- 
vantage of the mathematical simplifications associated 
with the existence of particle-hole symmetry. In addi- 
tion, one of the Lorentzian form is taken to add states 
to Nf,, while the other subtracts states so that there is 
conservation of total number of states when A''(e) is inte- 
grated over energy, i.e., AN{e)d€ — 0, where AN{e) 
is the change in the EDOS. The Lorentzian form has the 
important simplifying property that an energy integral 
in the definition of the conductivity can be done analyt- 
ically. 

For simplicity we also limit ourselves to the case of 
elastic impurity scattering. This case has been exten- 
sively studied in the approximation that N(e) is con- 
stant in the energy range about the Fermi energy which 
is significant for transport. For a constant Ni, the quasi- 
particle scattering rate l/rqp(ti;) = —Y,2{lu) is constant 
independent of energy uj. The conductivity takes on the 
well-known Drude form with constant transport scatter- 
ing rate which gives the half width of the Drude and is 
in fact equal to twice the quasiparticle scattering rate. 
When inelastic scattering is considered, the quasiparticle 
scattering rate can still be defined in terms of the self 
energy £2(0^), but noai, it acquires a temperature and 
frequency dependenceLHH. In this case, the quasiparticle 
scattering time rqp(w) is no longer equal to the optical 
scattering time rop(w), which is formally defined in terms 
of cr(w) — ai{Lu) + ia2{oj) through the formula 

Top(a;) 47r a{uj) Air al{uj) + al{uj)' 

where VLp is the plasma frequency which is related to the 
real part of the conductivity (Ji{uj) through the optical 
oscillator strength sum rule 

dcua.iu;)^^. (2) 

In contrast to the constant EDOS case, when N{e) 
varies with e around the Fermi energy, l/xqp and l/rop 
are no longer constant just as in the inelastic case and 
are not equal. Each acquires a separate dependence on 



energy. The imaginary part of the electron self energy 
I]2(w) becomes proportional to the self-consistent quasi- 
particle density of states N{u!) of the impure system. Im- 
purities broaden the pure crystal EDOS N{e) leading to 
N{uj). The optical scattering time defined in Eq. (|l|) also 
acquires oj dependence, and can ha quite different from 
the quasiparticle scattering timctHO, which is measured 
in ARPES experiments. Thus, optical and ARPES data 
give complementary information on What is mea- 

sured in ARPES is the single particle spectral density for 
a particular momentum k as a function of w. It is de- 
noted by A(k,uj) and is related to the self energy 
with k dependence suppressed, by 

The interpretation of optical results is now no longer 
straightforward. As an example of the complications 
that arise we note that for the case (as we will assume 
in most of our calculations here) when A^(e) conserves 
states, the integral over energy of the ARPES rate will 
also remain unchanged. This is because 1/Tqp(a;) is pro- 
portional to N{uj) and AN{e)de = is guaranteed 
when AN{e)de = 0. This sum rule however does not 
hold for l/rop(w), as have been previously discussednZl for 
the case of inelastic scattering processes and for the onset 
of superconductivity, with constant EDOS. Nor is the to- 
tal optical oscillator strength defined in Eq. @ constant. 
This arises because the integral over <7i(w) defining flp 
depends on an average of the EDOS N{e) around the 
Fermi energy over an energy scale defined by the impu- 
rity scattering and is not just dependent on A^(0) as it 
would be in the familiar constant EDOS case. There is 
also an attendant temperature dependence of Qp. As the 
temperature becomes comparable to the energy scale on 
which A^(e) varies significantly, the energy dependence in 
N{e) is effectively smeared out and we recover a simple 
Drude form. 

The energy dependent EDOS enters the formula for 
the conductivity in two places. First, the total current is 
the sum of the partial currents contributed by each state 
|k) in the electron system. When this sum is changed 
into an integral over energy a first factor of A'^(e) en- 
ters. But there is a second factor of A''(e) that also comes 
in from the quasiparticle scattering rate. This rate is 
proportional to the matrix element of the impurity po- 
tential which is to be averaged over all final states in 
which the electron can scatter. We can call this a final 
state effect. This second factor enters the ARPES rate 
which becomes proportional to the self-consistent N{uj). 
Clearly, ARPES and optical rates can no longer simply 
be proportional to each other. Both the initial and final 
states factors modify the optical scattering rate. 

In our calculations, we find that the factor of N{e) 
coming from the sum over partial currents from each 
electron has less of an effect on the energy dependence 
of <7i{uj) than does the modification of the underlying 
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ARPES rate due to final states effects. For a model of 
N{e) which has a depression in the EDOS at e = ep, 
which is, of course, compensated for at higher energies 
so as to conserve the total number of states, the first fac- 
tor of N(e) decreases the dc conductivity more than at 
finite frequency so that the overall effect is to lead to an 
apparent broadening of the Drude-like form for ai{uj). 
On the other hand, the ARPES rate is effectively re- 
duced at small cu by the final state factor of N{e). This 
sharpens the Drude-like line at small cu. Thus the two 
effects have opposite tendencies, compete against each 
other and partially cancel. In the specific cases consid- 
ered, the modifications in ai{uj) brought about by the 
changes in the ARPES rate are more important. 

In a final set of calculations we also consider the case of 
a step- function EDOS model. In a metal N{e) is expected 
to be finite at the Fermi energy although it could be small 
as compared to its value away from e — ep. With our 
step model we show that a small but finite value h of N{e) 
for |e| < some energy Eg about e = always leads to the 
existence of Drude-like peak in the optical response, in 
sharp contrast to the case h = when a gap forms and 
the Drude peak is completely eliminated. 

The paper is organized as follows. In Sec. II, we 
present a general theory of the optical conductivity in the 
case of impurity scattering. The simplest Drude limit is 
discussed in Sec. Ill as a reference. Sec. IV is devoted to 
a discussion of the effect of energy-dependent EDOS on 
optical conductivity, within a toy model for the EDOS 
involving two Lorentzian forms. This is followed by a 
parallel discussion, in Sec. V, for another EDOS model, 
the step model, which allows us to contrast metallic-like 
and semiconducting-like behavior. Finally, Sec. VI con- 
tains our conclusion. Some mathematics is shown in Ap- 
pendix A and B. 



II. GENERAL THEORY FOR OPTICAL 
CONDUCTIVITY 

In linear response theory, the optical conductivity can 
be expressed as 



(j{uj) = -[n{uj + id)~n{id)], 



(4) 



where Il{u!+i6) is the retarded polarization function. Un- 
der the assumpUan that vertex corrections are negligible, 
n(cc) + iS) read£3 

q2 poo ]\r(f\ r°° 
U{u; + ^5) = ^ dxf{x)A{k, x) 

x[G{(.,x + 00 + i5) + G[e,x - LO - id)], (5) 

where N{e)/Nb is the normalized EDOS with the 
constant background EDOS, f{x) the Fermi distribution 
function, G(e, x ± i5) the quasiparticle Green function, 
and f2po the bare plasma frequency which, for an energy- 
dependent EDOS, will be shown below to be different 



from the real plasma frequency Q,p defined in Eq. (|^). 
A(k, w) is the spectral density defined in Eq. (|3|). It is 
important at this point to emphasize that although, as 
we have stated, we have neglected corrections to the elec- 
tromagnetic vertex, the bare vertex itself can introduce 
further comphcations in Eq. (||). The EDOS factor iV(e) 
appearing in this equation comes from the conversion of 
a sum over all electron momenta into an integration over 
energy. But there is also a factor of the square of the 
Fermi velocity which is the electromagnetic vertex in our 
work, and this factor can have energy dependence. As 
we will not evaluate N{e) or for that matter the electron 
velocity v{e) from first principles but rather simply use 
a Lorentzian model, we can think that our model for the 
EDOS already contains the Fermi velocity and any de- 
pendence it may have on energy e. There is one caution 
we should make however. Later we will see, in our dis- 
cussion of the quasiparticle scattering rate, that a second 
factor of the EDOS enters and this one is not multiplied 
by the electron velocity squared. This second factor will 
further get renormalized by the impurity scattering, and 
so is replaced everywhere by the dressed EDOS N{lo). 
This should allow us to distinguish between these two 
factors of EDOS, and we will not emphasize this compli- 
cation further but it should be kept in mind. 

In the case of elastic scattering with no momentum 
dependence (no anisotropy), G(e,a;) can be written as 

G(e, a; ± «(5) — Gi(e, a;) =F i7rA(e, x) 

= [x - e - Y.i{x) ± ig{x)\ \ (6) 



where 



Gi(e,a;) 



X — Tii{x) — e 



[x-^,ix)-e]^ + g^x)-- 



(7) 



and g{x) = |E2(a;)| = Tqp^(a;) the quasiparticle scatter- 
ing rate. Si (a;) and g(x) satisfy the following Kramcrs- 
Kronig (KK) relation 



Si(a;) 



dx 



IT J —QQ X X 

Inserting Eq. (^) into Eq. (Q) leads to 



(8) 



47r J^^ uj 



X [A{€, x + Lo)- A{e, x~io)], 
^2(^) = ^/ I —f{x)Aie,x) 



(9) 



oo -^^> ^ — oo 



X [Gi(e, x + Lo) + Gi{e,x-io)-2 Gi(e, x)] ,(10) 

Eqs. (^ and ( p^ ) show how the two effects of the energy- 
dependent EDOS mentioned in the Introduction enter: 
one is the factor iV(e) coming from the sum over partial 
currents from each electron in the Fermi sea, the other 
arises from the electron spectral density factor A{'k,oo), 
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which contains a factor of the final states the particles 
are scattered into. 

The real and imaginary part of the conductivity, cri 
and (72, respectively, obey the KK relation, namely, 

(72 H = UJ / duj'a,{uj')P— -. (11) 

Eqs. ( ]ll| ) and (j^) lead to the useful relationship, 

(12) 



1 

lim <j2{uj) = — — ^. 
ixi— >oo 47r 



At zero T, Eq. (||) simplifies greatly and becomes 

N{e) f° dx 



02 

ai(Uj) = — TT , 



6 J-^ ^ 

Xv4(e, a;) yl(e, a; + uj). 
The dc conductivity immediately reads 



(13) 



a(0) 



TT 



47r 

1 

47r TT 



02 



A- 



dt 



N{e) 



(14) 



In Eq. (yl) we have used l]i(0) = 0. 

We assume that the impurity potential V is small and 
thus the impurity scattering can be treated within the 
Born approximation. Within the non-selfconsistent (nsc) 
Born approximation, the self energy reads 



■^Insc ^9nsc 

N{e) 



70 



de- 



(15) 



where 70 = riiV^Nb with rii the impurity density. 

For an energy-dependent EDOS, the selfconsistent (sc) 
Born approximation gives instead 



^iTHuj) = Ei,,(c^)-zg,,(c.) 

Nie) 



70 



de- 



G(e,tj + iO+), (16) 



where the full selfconsistent G appears on the right-hand 
side of Eq. ([l^), so this equation must be solved by 
successive iteration until convergence is achieved. 

Before leaving this section we emphasize that the fac- 
tor N{e) that enters Eqs. (^ and ^ is the EDOS 
when no impurities are present and comes from conver- 
sion of a sum over momentum to one over energy. It does 
not include any additional electron velocity factor v{e). 
This is to be contrasted with the formula for the conduc- 
tivity (Eqs. (|) and §)), in which we have suppressed 
a factor of ^(e)^, and so the EDOS factor that enters 
that formula should further contain any energy depen- 
dence there might be in ^(e)^. In our discussion below, 
we will describe the effect on the conductivity of each of 
these two factors separately. Note that gsc{i^) hi Eq. ( p^ ) 
is directly proportional to the renormalized EDOS N{uj) 
given by N{u;) = - / de[N{e)/Nb]A{e,Lu). 



III. DRUDE FORMULA FOR CONSTANT EDOS 

In the case of a constant EDOS N{e) = Nb, both the 
non-selfconsistent and selfconsistent Born approximation 
lead to Si = 0, g{uj) = tt^o = r/2. Eqs. (|) and (|l^) 
thus result in the well-known Drude formula 
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(Drudc)/ X "pO 



1/r 



(Drudc). N ^pO Uj/V^ 



47r i + {uj/ry 



(17) 
(18) 



^(Drudc)^^-j ( p^ ) is a Lorentzian function of uj with 

the half-width T. When Eqs. (|l^ and (|l|) are substi- 
tuted into Eq. (nl) which defines the optical scattering 
rate, we get t^^^ (uiY^'""^'''^ = 2g{uj) = T. In this simple 
case, the optical scattering rate is just equal to twice the 
quasiparticle scattering rate. Thus optical experiments 
access directly the microscopic information on the imag- 
inary part of the self energy. Further, Eq. (^7|) also gives 
the dc conductivity and the plasma frequency as 
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^(Drudc)cQN ^ ^ _ 

^ ' An T 



■pO- 



(19) 
(20) 



^(Drudc) j-^-j ^^ij-^-j (Drudc) as functions of are shown 
as green dot-dashed lines in Figs. ||, ||, || below and serve 
as a reference when we discuss the effects on the optical 
conductivity of an energy dependence in the EDOS. 



IV. OPTICAL CONDUCTIVITY FOR EDOS 
WITH TWO LORENTZIAN FORMS 

Now we are in the position to study the influence of 
the energy dependence of EDOS on the optical conduc- 
tivity. In this section, we consider the following model 
for EDOS, 



Nje) 
Nb 



1 



ai 



0.2 



n\ai + 



(21) 



where s > 0. The two Lorentzian forms of Eq. ( pi| ) guar- 
antee conservation of the total states: AN{e)de = 
J^^[N{e) — Nb]de = 0. For ai > 02, there is a hole, 
i.e., depletion of states, around the Fermi surface. This 
is shown as the black solid line in Fig. [|. While ai < 02 
corresponds to a peak, namely, additional states at and 
around the Fermi surface, as shown by the red dashed line 
in the same figure. The excess (missing) states are com- 
pensated for by a decrease (increase) in N{e) at higher 
energies beyond e ~ \fa\a2 ~ 7 in units of 70 for a\ — 5, 
di = 10 (oi = 10, a-i ~ 5]. 

From Eqs. (|l^) and ( |2l| ) we get the non-selfconsistent 
self energy based on the first-order Born approximation 



Slnsc('^) = 70 S 
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FIG. 1: Normalized density of states and ARPES scattering 
rate as functions of frequency in the Lorentzian-EDOS model 
shown in Eq. We used s = Stt. All energies are in units 

of 70. 



5nsc(w) = 7r7o 



7r7o 



TT 



ai 



0-2 



(22) 



so the black solid line and the red dashed line shown in 
Fig. |l| also represent gnsc/T^Jo vs cu. Eqs. ( p^ and ( |2l| ) 
yield the two equations determining the selfconsistent self 
energy 



Eisc(t^) = 70S 



LUq 



LOq 



all 



9I2 



gsc(w) = 7170 + 70S 



9al 



9a2 



9I1 



^0 + 9I2 



(23) 
, (24) 



where = uj - Sisc(^), 9ai = 3sc('^) + «!, and ga2 = 
(7sc(w) + 02- Eqs. ( p3| ) and (|2j) can be easily solved nu- 
merically. There exists an extensive literature on the 
effect of energy dependence of the EElQS 
tron self energy and on other propertiesElrH 
Some aspects of the superconducting state as well as nor- 
mal state have been explored although much of the liter- 
ature deals with a peak, while here we have emphasized 
a pseudogap. The imaginary part Qsc obtained from Eqs. 
( p3| ) and ( p^ ) for ai > 02 and for a2 > ai are shown in 
Fig. as magenta long-dashed line and blue dot-dashed 
line, respectively. Comparing g^sc and g^c in Fig. |l|, we 
see that selfconsistency smoothes out the ARPES rate 
because it is broadened out by impurity scattering. In 
fact, this rate is simply proportional to the fully renor- 
mahzed EDOS N{e) defined by Eq. (U) with factor 7770 
left out. Eq. (|2^) has the same form as Eq. ( pT| ) with 
the broadening in the Lorentzian form, i.e., ai and 02 
replaced by ai + gsci^o) and 02 -I- gsc('j-'), respectively. An 



iterative solution is required until convergence is reached. 
In Fig. [I] we see that for the case of the hole (peak) at the 
Fermi surface N{e)/Ni, at e = is 0.5 (1.5) and the appli- 
cation of selfconsistency changes these numbers to about 
0.7 (1.2). The smearing due to the impurity scattering is 
considerable in the EDOS. 



A. Constant ARPES rate g{uj) = r/2 but 
energy-dependent N(e) 

As we mentioned in the Introduction, there are two 
effects of EDOS on the optical conductivity, one is from 
summing over partial currents. This shows up as the ex- 
plicit factor of N{e) in Eqs. ^ and (|l^), the other is 
the final-state effect entering in the ARPES rate which 
shows up in Eq. (|2j). To see clearly the different roles 
these two effects play, we first switch off the effect of 
energy-dependent ARPES rate, and replace it by a con- 
stant g{uj) = r/2. Correspondingly, Si = 0. 

We focus on the T = case, because it is simplest and 
allows us to produce partially analytic results. Eqs. ( p^ ) 
and ( p^ ) immediately yield 



0-1 (cj) 

0-2 (w) 



An TT 
An n 



where 



m ^1^/2 U/2' r/2 j'^^^^ 
de ^ 7('-) f — (26) 



1 



o _ , , o , ^ 1 arctan(u) + e) 
72 w(w2 _|_ 4) 1 

+ e)2 + 1 



UJ 



1 



1 

72 (Ii([I)2 + 4) 

4 



In 



{Cu + ef + 1 



e2 -f 1 

— arctan(tl' + e) H — ^ 

a; e- + 1 



(27) 



(28) 



The expressions for the dc conductivity and plasma fre- 
quency become particularly simple and very revealing. 
From Eqs. (Ilj) and (||) we get 



a(0) 



02 



2 

An f 



02 



where 



n (a;2 



^ A^(xr/2) 

dx Yo{x), 

'^.^r.(.), 



Y,{x) = -- 



1 



n x^ 



1 



(29) 
(30) 

(31) 



We have plotted the two functions Yq{x) and Yi(x) in 
Fig. H. Both peak at a; = 0, and decay rapidly on a scale 
of e equal to a few times r/2. For small F, both (t(0) 
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and f2p depend strongly on the EDOS at and around 
the Fermi surface. If N{e ~ 0) < Nh 



for ai > a2, cr(0) < cr(°™'^'=)(0) and n„ < n„n immedi- 



which is the case 

pO 



ately follow. In the opposite case for ai < a2 in which 
N{e ~ 0) > Nb, we get cr(0) > cr^D^"'^'^) (0) and flp > flpQ. 
Besides, the peak in Yi{x) is broader than that in Yo{x), 
implying that the dc conductivity and the plasma fre- 
quency cr(0) and f2p do not scale with each other for 
energy-dependent EDOS. The region in energy around 
the Fermi energy that is most important in determining 
cr(0) and fi^ is given by F. If, as we have assumed so far, 
the scale for F is much less than the energy scale that 
controls important variations in the EDOS which in our 
case is (01,02) in the Lorentzian form, then it is mainly 
the value of N{e) at e = which comes in. But when F 
is of the same order as (ai, 02), this is no longer the case 
and the details of the variations in N{e) are importantly 
sampled. Finally when F is much greater than (oi, 02), it 
will be only the size of the background Ni, that matters. 
This is the limit in which we regain the simple Drude 
model of Eqs. @ and 

It is easy to find from Eqs. ( ^ ) and ( po| ) that 
Imiui^nc o'2(w) = ^p/47rcj, indicating the KK relation Eq. 



UJ) is automatically obeyed as we expect. For a general 
value of w, (Ti and 172 in Eqs. ( p5[ ) and (|2^) need to be 
obtained numerically. Inserting results of Eqs. (^5|), ( |2^ ) 
and (^0|) into Eq. (|l|) allows us to obtain the optical 
scattering rate. Results will be discussed in Sec. 3. C. 



where D = D/g{x), a;o = [a: - T,i{x)]/g{x), g± = g{x ± 
^)/9{x): x± = [x±w-Ei(a;±w)]/g(x), and 



and J^g'™-* defined in Eqs. (AI- A3[) in Appendix A. 
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FIG. 2: Scaling functions Yo{x) = 2/[ti{x^ + if] and Yx{x) 
2/[i.{x^ + l)]. 



C. Energy-dependent N{e) and ARPES rate 



B. Constant A''(e) but energy-dependent ARPES 
rate g{ui) 

If instead we switch off the energy-dependent EDOS 
by replacing it with N{e) — Nj,, but turn on the energy 
dependent ARPES rate g{uj) alone, we are able to see the 
second effect on a, coming from the ARPES rate alone. 
We have referred to this aspect of the problem earlier 
as a final state effect. Since in the present work the self 
energy is momentum independent, it is convenient to first 
integrate over e in Eqs. (|l|) and (|l^) at T = 0. After 
some algebra we obtain 



(cons) / N 



^lo r° dx 1 



lini r ^J— 

in TT D^oaJ_^ LO g(x) 

xT^''\D,g+,xo,x+) 
.9+ + 1 



47r UJ g{x) {g+ + If + {xq - x+f 

(cons)/ X ^pO 1 f° dx 1 



(32) 



47r TT D^ooJ_^ UJ g{x) 

[ J2 ^^"''HD, g^, XO, X,) - J't\D, Xo)] 



^lo dx 1 



x+ - Xq 



47r UJ g{x) {g+ + ly + [xq - x+f 



(33) 



If we want to include in the calculations both the di- 
rect factor of N{e) in Eq. (^ and the energy-dependent 
ARPES rate so as to study the competition between the 
initial-state effect and the final-state effect, we need to 
return to Eqs. (|2l]) and (|2^) or (p5) and insert these into 



Eqs. ( [l3|) and to find (Ji{uj) and t~^{uj). Note that 
the Lorentzian forms used for N{e) allow us to do the 
integral over e analytically, which greatly simplifies the 
numerical work. The expressions for (Ji{uj) and 172 (w) af- 
ter integrating over e are given in Appendix B. They are 
complicated but can easily be handled numerically. 

Our first numerical results are presented in Fig. H, 
which has two complementary frames (a) and (b). In 
all curves oi — 10, 02 = 5 with s = Stt, in units of 
70- This corresponds to a depression in the EDOS (See 
Eq. (^)) over its background value, as drawn in Fig. |l| 
solid black curve with A^(0) reduced by a factor of two 
compared to the background value. The missing states 
at small e show up at larger e and there is conservation 
of states by arrangement. Although electronic states are 
preserved this does not mean that the corresponding op- 
tical oscillator strength is, as we will see. In Fig. |^ (a) 
we show five curves. All give the real part of the optical 
conductivity (Ji(uj) in units of f2po/47r as a function of 
energy uj in units of 70. The first in green (dot-dashed) 
is for reference and is the usual Drude form (given in Eq. 
(p7|)) that obtains when the EDOS is constant as well as 
the elastic scattering rate. In this case, the width of the 
Drude is simply 27r7o = F and the oscillator sum rule 
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FIG. 3: (a) Real part of optical conductivity and (b) Optical 
scattering rate as functions of frequency in the Lorentzian- 
EDOS model shown in Eq. ( ^l| ) with ai — 10, 02 = 5 and 
s = 5n. We used 70 = 1. Labels of curves in (b) mean the 
same as in (a). 



gives r2po/8. This no longer holds when energy depen- 
dence is included in the EDOS. Introduction of such a 
dependence in the theory leads to several modifications 
and we will take these in steps of added complications. 
The simplest modification is that an overall EDOS factor 
enters the sum over partial currents involving each par- 
ticipating electron. This was referred to previously as the 
sum over initial states. Including only this factor with a 
constant approximation for the scattering rate g{i-j) — tt 
in units of 70 gives the blue curve (long-dashed). The 
reduction in N{£) at e = by a factor of 2 translates into 
a substantial reduction in aiico) at w = although by 
a factor that is substantially less than 2. While the line 
shape is no longer perfectly Drude, its width at half max- 
imum has increased over the Drude case T = 27r7o and 
could lead one to conclude that the optical scattering'Tate 



has increased. It is already clear from this remark that 
optical and quasiparticle scattering rate are no longer 
the same. In fact, this is shown explicitly in Fig. || (b). 
The green curve (dot-dashed) is constant but the blue 
curve (long-dashed) now exhibits a slight energy depen- 
dence and is everywhere larger than twice the quasiparti- 
cle rate. Another important modification brought about 
by the introduction of an energy dependence in N{e) is 
that the optical sum rule defined in Eq. (^) as the sum 
over the real part of the conductivity has the plasma 
frequency reduced from f2po to Qp = 0.875 QpQ. Note 
that this occurs although no states are lost in N{e). On 
the other hand, when the initial-state EDOS factor is 
taken to be independent of energy (constant Ni,) there 
is no change in plasma frequency even if the quasiparti- 
cle scattering rate is energy dependent as shown in the 
red solid thin curve which was computed for constant 
N{e) = Nh, but with Y,nsc of Eq. (^) obtained in a non- 
selfconsistent theory for the quasiparticle scattering rate. 
By non-selfconsistent we mean a first iteration of the self 
energy equation in Eq. (16) as defined in Eq. (|l^). This 
modulates the scattering rate with precisely the same 
EDOS factor N{e) that we used for the initial state sum 
in the blue curve. We see that in a real sense this factor 
has the opposite effect on the shape of the conductiv- 
ity <7i{lu) vs lu in that it increases its value at w = 0, 
beyond what it is in the pure Drude case, and also ef- 
fectively sharpens up the curve. This can be seen more 
quantitatively in Fig. ^ (b) where the red solid thin curve 
for the optical scattering rate falls everywhere below the 
green curve (dashed-dotted) of the Drude theory and also 
even further below the blue curve (long-dashed). Com- 
bining these two effects brings us back closer to our orig- 
inal Drude than including each separately, at low fre- 
quencies where the effects are biggest. This expectation 
is born out in the solid black curve which includes ini- 
tial state iV(e) factor and non-selfconsistent quasiparticle 
scattering rate (?nsc(^)- Note that the plasma frequency 
(rip/ripo — 0.848) is not changed much from its value 
of 0.875 in the blue curve, which shows that the plasma 
frequency is mainly sensitive to the initial N{e) and is 
less sensitive to the details of the quasiparticle scattering 
rate. The corresponding optical scattering rate is shown 
as the solid black curve in Fig. ^ (b) . This curve contrasts 
greatly with the previous two curves blue (long-dashed) 
and red (solid thin). In both these curves there is no com- 
pensation in the scattering rate as compared to the Drude 
case in the sense that the blue curve is always above and 
the red curve always below. By contrast, the black curve 
(solid) is below at small lu and above at larger oj. While 
there is some cancellation, no sum rule applies to the area 
under l/rop(ci^) when it is integrated over lu. This rep- 
resents a real difference between quasiparticle scattering 
rate and optical scattering rate since for the ARPES rate 
a sum rule does apply which is directly connected to the 
sum rule on N(e). This shows clearly that, while optics 
can give microscopic information on scattering rates, it is 
not easy to relate the information so obtained, with the 
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characteristics of the self energy which is, in the end, the 
fundamental quantity and is the quantity we would like 
to measure directly. 
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FIG. 4: (a) Real part of optical conductivity and (b) Optical 
scattering rate as functions of frequency in the Lorentzian- 
EDOS model shown in Eq. with ai = 5, 02 = 10 and 

s = 5n. We used 70 = 1. Labels of curves in (b) mean the 
same as in (a). 



One final element of importance has been neglected 
so far. To get the effective scattering rate in a system 
with energy dependent EDOS with impurities, it is nec- 
essary to solve the self energy equation ( |l6|) selfconsis- 
tently through repeated iteration of Eqs. (p3[) and (24). 
Impurities will smear out a valley in N{e) and change 



it to its selfconsistent value N{e), which is what finally 
enters in the impurity scattering. When this is done we 
obtain a final curve for cri(w) vs uj, the magenta curve 
(dashed) of Fig. ^ (a) and for the optical scattering rate 
in Fig. g (b). As we would have expected, the self consis- 
tency smoothes out the curves but this does not introduce 
any new physics. The value of the plasma frequency is 



also not changed much over its non-selfconsistent value. 

Before moving on to the case of a peak in the EDOS 
at the Fermi energy, we bring up once more the compli- 
cation that arises from the electromagnetic vertex which 
introduces a product of two electron velocity w(e)^ in the 
formula for the conductivity which effectively introduces 
a further energy dependence in the factor of N{e) ap- 
pearing directly in Eq. (0) . This factor is not present 
in the self energy Eqs. ( |l5| ) (non-selfconsistent) and ( |l^ ) 
(selfconsistent). This means that the magenta (dashed) 
curve of Fig. |3| (a) could be further modified through 
the introduction of a w(e)^ factor in our model N{e). If 
we look at band structure calculations, we note that the 
product of AL(e)?;(e)^ is often less dependent on energy 
than is A^(e)c3. This means that in this case, the final re- 
sults for the conductivity might move some way towards 
the results of a selfconsistent theory for the self energy 
with constant N{e) as the explicit factor of EDOS in Eq. 
(i). One curve shown in Fig. |^ (a) is for N{e) in Eq. 
(I) equal to a constant with non-selfconsistent self en- 
ergy (red dotted line). The selfconsistent case would be 
smoothed out a little more as compared with this curve. 
Of course, in this work, we do not wish to tie ourselves 
to band structure calculations but rather treat a general 
model for the energy dependence. 

In Fig. U we show results when there is a peak in N{e) 
rather than a depression or gap. No new physics is associ- 
ated with these curves. In a sense, the effects are opposite 
of those described in some detail for the previous case. 
For example, comparison of the results of full calcula- 
tions including the self energy computed selfconsistently, 
shows that cti {uj) is above the simple Drude curve re- 
flecting the fact that there is a peak in N{e) centered at 
e = 0, which enhances the value of (Ti(lu) around uj = 0, 
because there are more electronic states around the Fermi 
energy. Also when this same effect is included in the scat- 
tering process the scattering rate is increased at small u) 
as compared to its Drude value, which makes the curve 
ai{u)) vs to look broader. In the example considered, the 
actual frequency dependence of the selfconsistent optical 
scattering rate is not very large (see the magenta dashed 
curve in Fig. ^ (b)) and its variation is smaller than found 
for the self-consistent quasiparticle scattering rate shown 
as the blue dot-dashed line in Fig. |l]. These two curves 
are similar but are not the same in detail. In particular, 
we stress again that there is a sum rule which holds for 
the quasiparticle rate as compared to its Drude counter- 
part, i.e., no such sum rule applied to the area under the 
optical rate. 

Another consequence of the energy dependence in 
N(e) is that the conductivity will change with temper- 
atures, an effect which does not arise in the ordinary 
Drude case. This is illustrated in Fig. ^ where we show 
the real part of the conductivity cri(w) vs u for three 
cases. The black solid curve repeats our previous re- 
sults for the case of a gap at the Fermi energy and non- 
selfconsistent ARPES rate, shown also as the black solid 
line in Fig. ^ (a). The green dot-dashed curve is the ordi- 
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FIG. 5: Real part of optical conductivity at different temper- 
atures in the Lorentzian-EDOS model shown in Eq. . We 
used the same ai, a2, s and 70 as in Fig. H. 
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nary Drude for comparison. The red dashed curve shows 
how the solid curve evolves with increasing temperature. 
It represents results of (Ti (uj) vs lo with non-selfconsistent 
self energy for temperature T = 10 in units of 70 which 
is related to the quasiparticle scattering rate. It is clear 
that the evolution is towards restoring the curve to its 
simple Drude value as can be expected when tempera- 
ture or impurity effects are sufficiently strong that they 
wash out the energy dependence in A^(e). Note that for 
the case considered here, the energy scale for the struc- 
ture in N{e) is 10 in our units and this scale is of the same 
order as the temperature. Similar smearing effects are ex- 
pected when impurity scattering is increased sufficiently 
that the impurity scattering rate becomes comparable to 
the energy scale of the structure in A^(e). In this case we 
do not show a curve analogous to Fig. ^, but instead we 
show the change in the optical oscillator strength under 
the curve, i.e., the plasma frequency. Before presenting 
these results we point out that in Fig. || at T = the 
plasma frequency ilp/ilpo = 0.848, while at T = 10 it 
has moved up to 0.976 close to the simple Drude case, 
so that effects of energy dependence in N{e) are pretty 
well washed out in Qp as they are in the full ci (cj) vs uj 
curve. In Fig. ^ we present equivalent results for the im- 
purity scattering and compare with temperature. What 
is shown is the plasma frequency as a function of impu- 
rity scattering rate (or temperature). The red dashed 
and blue dot-dashed lines are for a peak in N{e), and the 
black solid and magenta long-dashed lines for a valley at 
the Fermi energy. We see that temperature and impurity 
have very similar effects on the plasma frequency and 
that when the scale on the horizontal axis is of the order 
of the scale defining the structure in iV(e), we recover in 
both cases the Drude plasma frequency as we expected. 

The above theoretical calculations apply to materials 
in which a pseudogap exists in the electronic band struc- 



ture. The experiments need to be done at relatively low 
temperatures where the electrons are dominantly scat- 
tered by impurities and the systems must remain in the 
normal state. The underdoped cuprates exhibit a pseu- 
dogap of size 300 K or so. This pseudogap has its origin 
in many-body effects as we mentioned in the Introduc- 
tion. Nevertheless, in as much as the effective EDOS 
around the Fermi surface is depressed, our results should 
give some qualitative guidance about the physical conse- 
quences of such a EDOS depression. Of course, in the 
cuprates a superconducting transition occurs, but Tc can 
be low, of order 20 K. In this intermediate temperature 
regime, however, the inelastic scattering due to electronic 
interaction may already be significant, which complicates 
the comparison. Disordering in the Cu02 planes, such as 
replacing copper sites by nickel, in the underdoped high- 
Tc cuprates, will efficiently enhance the elastic scattering, 
so that it becomes dominant. Alternatively, a magnetic 
field can be used to quench the superconductivity, en- 
abling measurements at much lower temperatures. 



V. OPTICAL CONDUCTIVITY IN "STEP 
MODEL" 

Another EDOS model which allows us to get simple 
analytical results and also helps in developing insight into 
the effect of N{e) on (j{uj) is the so-called step model. 



h, 
1, 



|e| < Eg 
kl > Ea 



The non-selfconsistent self energy from Eq. 

-E„ 



^Insc 
ffnsc(w) 
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The selfconsistent self energy from Eq. (|l^) is determined 
from the foUowing equations 



Sisc(w) = (1 - /i) Y In 



(wo + ^3)2+.g2,M 



gsc{'^) = 7o ^ ^ (1 ^ ^) arctan 

f UJO~ Eg 

— arctan — — ^ 

V 5sc(w) 



^0 + -Eg 

.9sc(^^) 



(37) 



(38) 



where loq — ~ Sisc(^)- Eqs- ( |37| ) and ( |3S| ) need to be 
solved numerically. 

To calculate cri{uj) in Eq. j]!^), we can first integrate 
over e. Inserting Eq. (|3^) into Eq. (|l3|), we find that 



(71 (w) 



1 /■° dx 
47r TT 



^(.9 + 1) 



uj { g{x) {g + 1)2 + [xo - x)^ 



-(l-/.).F("'(£;g,5,xo, 



(39) 

7(a; + w)/g(a;), 
-w)]/g(x), and 
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where is found in Appendix A, ^ = 

xq = [x — T,i(x)]l g{x), X = [x + ijj — Y,i(x- 
Eg = Eg/g{x). 

In Fig. R we show the numerical results that we have 
obtained in the step model. Fig. |7| (a) gives the quasi- 
particle scattering rate g{LLj) vs cu for four different cases. 
The green dot-dashed line is the non-selfconsistent result 
for g{uj) in the semiconductor-like model, i.e., h = 0. As 
indicated in Eq. (^6|), g{io) is exactly zero until w = 10 
and jumps to tt for a; > 10. This curve is for compar- 
ison with the full selfconsistent quasiparticle scattering 
rate obtained from iterating to convergence Eqs 
and (|8|). The numerical results are the blue long-das 
curve. We see a sharp scattering edge at flo — 6.1 below 
which the scattering rate is zero, so that now the sharp 
onset has moved to lower energy, as compared to the 
non-selfconsistent case. Note that g(w) is proportional to 
the fully renormalized (from impurity scattering) EDOS 
N{uj) as defined before and this shows smearing of the 
step function edge that existed in the non-selfconsistent 
case. The impurity scattering also accounts for the reduc- 
tion in scattering above to — 10. In fact, there is a sum 
rule on N{uj) so that the area under each curve, selfcon- 
sistent and non-selfconsistent, remains unchanged. So, 
the increase in fully renormalized EDOS below w = 10 
is fully compensated for at larger w. This also holds 
for the second set of two curves. In this case h — 0.1 
instead of being zero. The red dashed curve gives the 
non-selfconsistent results while the black solid curve is 
for the selfconsistent case. In this case, there is always 
a finite EDOS at all frequencies, but it still rises sharply 
at w = 10 in the non-selfconsistent case and at w ~ 6.1 
in the selfconsistent calculations. However, as compared 
with the h = case, the rise is not sharp. There is some 
rounding of the edge. Conservation of states applies how- 
ever. 

The corresponding results for the real part of the con- 
ductivity ai{uj) vs UJ are shown in Fig. |^ (b). We see, as 
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FIG. 7: (a) ARPES rate and (b) Real part of the optical con- 
ductivity. as functions of frequency in the step model shown 
in Eq. (p4). We used Eg = 10 and 70 = 1. 



we expected, that for the pure semiconductor-like model, 
there is no optical absorption until uj = 2Eg = 20 in the 
non-selfconsistent case (green dot-dashed curve). This 
also holds for the selfconsistent case (blue long-dashed 
curve), but now the edge has moved to oj = 2^0 ~ 12.2. 
The other two curves are metallic in nature because now 
h — 0.03 7^ 0. The finite value of h at small w guarantees 
that the response in this region is metallic-like and we 
see a sharp Drude-like peak centered around cu — 0. At 
higher energy the semiconductor-like behavior of the two 
previous curves remains although the main rise in the 
conductivity has been shifted to lower uj as compared 
with u — 20 and 12.2 respectively for non-selfconsistent 
and selfconsistent case with h = 0. Between the Drude 
response at small lo and semiconducting-like response at 
high UJ the curve gets filled in and shows sharp structures 
corresponding to the sharp step assumed to exist in the 
initial EDOS N{e). 
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We next contrast more sharply the two qualita- 
tively different limits discussed above using an analytical 
method to get simple although limited results. 

Case 1. h — 0: Semiconducting-like behavior. 

In this case, g{uj) ~ for w < Op. (Note that for the 
non-selfconsistent ARPES rate, fig = Eg, and for the 
selfconsistent ARPES rate, fio < E„.) 

Now we look at cri(ijj) in Eq. ([l3|). Noting that —uj < 
x < 0, we obtain: 

i) For < Lu < Qo, both A(e, x) and A{e, x + uj) are 
zero, leading to (Ti(w) = 0; 

ii) For fio < < 2f2o, A{x + uj) = for -uj < x < 
— rio, while A{x) = for — fio < x < 0. Therefore, 
ai{uj) = 0. 

To conclude, <Ti{uj) = for uj < 2ilo. 

Case 2. h — 0+: Metallic-like behavior. 

In this case, we expect to see the Drude peak at small 
UJ. This can be nicely shown in the dc conductivity. By 
substituting Eq. (34) into Eq. ( [T^ ) we obtain 
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It is clear that iax h — \, Eq. ( [lO[ ) recovers the Drude 
result, CT^i^^^'^HO) = (f&)/47r)/2g(0). In the limit oih^ 
0+ , we get from Eq. (|o|) that 
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For the non-selfconsistent case, (7nsc(0) — T^"foh, thus 
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(42) 



For the selfconsistent case, however, Eq. (|3^) produces, 
in hnear order of h, gsc{0) — T^loh/[l — 27o/£'g) (for 
270 < Eg), leading to 
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1 



An 2n-fo 



1-1^ I, (43, 



which is different from the corresponding Drude value. 

We can also compare the above results obtained in 
the EDOS-step model for a finite h value with those ob- 
tained in the EDOS model with one Lorentzian form, 
N{e)/Nb = 1 - {s/n)a2/ial + e^). We find that in the 
Lorentzian- form model, ai{uj) shows a sharp Drude-like 
peak at low frequencies and semiconducting-like rise at 
around u a2, which are similar to what is shown in 
the solid black curve in Fig. ^ However, since both the 
EDOS and the ARPES rate are smooth in this model, 
the conductivity curve is always smoothly evolving as a 
function of frequency, and the semiconducting-like rise 
can never be as sharp as in the step model. 



VI. CONCLUSION 

We have studied the effect, energy dependence in the 
EDOS around the Fermi energy, has on the frequency de- 
pendent conductivity, and the derived optical scattering 
rate. Only elastic impurity scattering is accounted for. 
To keep the mathematics simple, we employ for N{e) a 
constant background modified by two Lorentzians both 
centered at the Fermi energy so as to provide particle- 
hole symmetry. One adds states the other subtracts, so 
that the sum over all states is not modified. One of the in- 
tegrals defining the conductivity in terms of the electron 
self energy S(w) can be done analytically. This greatly 
simplifies the numerical work and allows us to find cer- 
tain important limits in closed form. The conductivity 
is modified in two important ways. One comes from the 
energy-dependent EDOS factor which enters when a sum 
is to be carried out over all the electrons that contribute 
to the current. We referred to this as the initial state 
factor. The second factor comes from modification of 
the quasiparticle scattering rate which is no longer con- 
stant and in fact becomes proportional to the selfcon- 
sistent EDOS N{e). This quantity differs from N{e) in 
that it accounts for the smearing of the EDOS brought 
about by the impurities. It is to be computed from the 
selfconsistent self energy equations. Both factors have a 
profound effect on the conductivity, but their individual 
contribution cannot easily be separated. For a depres- 
sion in N{e) at the Fermi energy, the first factor reduces 
the real part of the conductivity at and around uj = 
(dc conductivity) effectively, making the resulting opti- 
cal scattering rate appear larger than its value in the pure 
Drude model, i.e., the width of the curve at half max- 
imum is increased. On the other hand, the second fac- 
tor acts in the opposite manner. Since the quasiparticle 
scattering rate becomes proportional to N{uj), this quan- 
tity now acquires a frequency dependence which, in the 
model under consideration, reduces the scattering rate 
at and around uj = 0, because N{uj) is less than the 
constant background value there. This has the effect of 
increasing the height of the Drude at small uj and mak- 
ing it narrower as compared to the constant case. The 
net result of both effects in the selfconsistent case is to 
produce a curve for cri(w) vs uj which is reduced below 
the pure Drude at small uj. The corresponding optical 
scattering rate is depressed at small uj below its constant 
Drude value and then becomes larger at higher uj. This 
mimics the variation in the quasiparticle scattering rate 
which is proportional to N{uj), but there are important 
differences. In particular, while a sum rule applies to the 
ARPES scattering rate which reflects quite directly the 
sum rule on the EDOS, the optical scattering rate dis- 
plays no such property. It is clear then that when there is 
important energy dependence in the EDOS at the Fermi 
energy, optical and quasiparticle scattering rates are not 
as simply related as in the Drude model. This compli- 
cates the process of extracting microscopic information 
from optical scattering rate data. What one ultimately 
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wants is detail information on the self energy. 

There are several other complications that arise that 
need to be commented upon. The plasma frequency 
which gives the optical oscillator strength, i.e., the area 
under the real part of the conductivity, becomes depen- 
dent on temperature and on impurity concentration. We 
show that the plasma frequency depends on a range of 
states around the Fermi energy with the energy scale as- 
sociated with this range given by the impurity scatter- 
ing rate. It is not just its value at the Fermi surface 
which matters. As the impurity rate is changed, the 
range of important values of N{e) is also changed and 
so is i7p. In the models considered this can be a signif- 
icant effect. Temperature can also smear out the region 
in energy around the Fermi surface and so also impact 
on the value of U,p. As temperature is increased, we find 
that the EDOS effects become gradually less important 
and the entire curve for ai{Lj) vs u) moves towards the 
Drude form with constant background only in N{e). An- 
other result of interest is that the region in energy most 
important in determining the plasma frequency is differ- 
ent from that determining the dc conductivity, so that 
these two quantities will not respond in the same way 
with change in impurity content. 

Finally we have considered a step function model for 
N{e) around the Fermi surface with N{e) — hNj, for < 
e < Eg, and N{e) = Nt for e > Eg, where Eg is a gap 



energy. For h = 0, this is a semiconductor-like model 
which leads directly to zero conductivity in the range < 
u! < 2ilo (ilo — Eg for non-selfconsistent ARPES rate, 
and flo < Eg for selfconsistent ARPES rate). We find 
that for finite h, however small it may be, the situation 
is radically different and an intrinsically metallic behavior 
is always obtained. At small uj, there is a very narrow 
Drude-like peak followed by a depression, and at higher 
energies u; > ^Iq a semiconductor-like behavior is again 
observed. The value of the dc conductivity is unchanged 
from its value with N{e) = N^, everywhere provided h 
1. Here we have described only the case when the self 
energy is treated in a first iteration. We find however 
that the situation shows no qualitative change when a self 
consistent theory is considered except for the important 
difference that the gap flo is reduced by the interactions 
and the dc conductivity has a correction linear in jo/Eg. 
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APPENDIX A: EXPRESSIONS FOR Jf("^'=>, JF^'"^' AND T^'""^ IN EQS. @ AND 



As mentioned in Sec. III. B, in the case of a constant N{e) = Nb and energy-dependent ARPES rate, we can first 
integrate over e. The resulting optical conductivity at T = is shown in Eqs. (B2h and (p3), where 



^(-)(^,5,xo,x) =do [5(2/'+.9'-l)/i + (2/'-.9' + l)/2-5 2//3], (Al) 
g, xo,x) = ^[- 2y(y' + + l)/i + 4 g y f2 + {1 - g' + y^)h] , (A2) 
J't'\D, xo) = 4x0^ [{xo - D f + 1] [(xo + bf + l] "\ (A3) 

with do — [(g -I- 1)^ -I- y^] ^ [(g — 1)^ -I- y^] ^ , y = xq — x, fi = arctan(i) + xq) + arctan(l) — xq), /2 = arctan [{D + 
x)/g] + arctan [{D - x)/g] , and /a = ln{ [{D -xof + l]/ [{D + xo)^ + l] } + ln{ [{D + xf + 5^] / [(^ - xf +52]}. 



APPENDIX B: OPTICAL CONDUCTIVITY FOR ENERGY-DEPENDENT ARPES RATE IN THE 

LORENTZIAN DOS MODEL 



For the EDOS model with Lorcntzian forms shown in Eq. (21) and the energy-dependent ARPES rate, we can 
first integrate e to simplify the numerical work. The resulting conductivity reads 



ai{uj) = CTj™"*'''(tj) -I- 5ai{uj,ai) - ^1(^,02), 

6a2 (^^,02), 



(Bl) 
(B2) 



n (cons) 1 (cons) i . 

where a\ and a 2 are shown m Eqs. 



5iji{lo, a) 



47r 27r 



dx 

UJ 



and (|33|), and 

Xo Xi + ga+gia- 



Xo X2 - ga+92a+ 



{xl + gl+) {xi + <?L-) {4 + 9I+) {A + 9laA 
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_^^pa so. dx 



2{xi - xo) xigi + (.go + (a^ + xf ~ gf) 



47r TT w [{xi - xo)^ + {go + 51)2] [{a^ 



glr+4xlgiy 



(B3) 



. , . s r° dx 

47r 2n J _^ w 



a;26'a- 



+ 



+ 



xogia+ + Xiga+ 



{xl + gl_) {xl + 9^+) [xl + gl+) {xf + gl^_) {xj + 5^+) (x? + g^^) 



+ 



Xog2a+ + X2ga+ "^Xg^-i 

{xl + gl+){xl + 9L+) {x^ 



sa dx 

47r TT w 



(xq - X2){a? +xl- gl) - 2xogo{go + ^2) 

[(a:o - + {go + 52)2] [(^2 + ^2 _ ^2)2 + 4^2^2] 



(xi - .To)(a^ +xl ~ gl) - 2xigi{go + gi) 



[(xi-xo)2 + (50+51)2] [{a^ 
with Xo = X — Si(x), Xl = X + w — Ei(x + w), X2 

5a± =50 + 0, 5la± = 51 ± a, 52a± = 02 ± a. 



(B4) 



c - w - Si(x - w), fifQ = g{x), gi = g{x + u)), 52 = g{x - w), 
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